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Open 



A substantial part of the Earths' soil organic carbon (SOC) is stored in Arctic permafrost peatlands, 
which represent large potential sources for increased emissions of the greenhouse gases CH 4 and 
C0 2 in a warming climate. The microbial communities and their genetic repertoire involved in the 
breakdown and mineralisation of SOC in these soils are, however, poorly understood. In this study, 
we applied a combined metagenomic and metatranscriptomic approach on two Arctic peat soils to 
investigate the identity and the gene pool of the microbiota driving the SOC degradation in the 
seasonally thawed active layers. A large and diverse set of genes encoding plant polymer-degrading 
enzymes was found, comparable to microbiotas from temperate and subtropical soils. This indicates 
that the metabolic potential for SOC degradation in Arctic peat is not different from that of other 
climatic zones. The majority of these genes were assigned to three bacterial phyla, Actinobacteria, 
Verrucomicrobia and Bacteroidetes. Anaerobic metabolic pathways and the fraction of methano- 
genic archaea increased with peat depth, evident for a gradual transition from aerobic to anaerobic 
lifestyles. A population of CH 4 -oxidising bacteria closely related to Methylobacter tundripaludum 
was the dominating active group of methanotrophs. Based on the in-depth characterisation of the 
microbes and their genes, we conclude that these Arctic peat soils will turn into C0 2 sources owing 
to increased active layer depth and prolonged growing season. However, the extent of future CH 4 
emissions will critically depend on the response of the methanotrophic bacteria. 
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Introduction 

Recent estimates point out that 277 Pg of soil organic 
carbon (SOC) are stored in Arctic peatlands 
(Tarnocai et ah, 2009), which corresponds to 1/3 of 
the C0 2 in the atmosphere. These peat soils have 
acted as carbon sinks since the Holocene (Post et ah, 
1982; Harden et ah, 1992). In contrast, they are 
substantial sources of methane (CH 4 ), releasing 
~35Tg per year, ~6% of the global CH 4 emissions 
(Cao et ah, 1996). CH 4 is a more potent greenhouse 
gas than carbon dioxide (C0 2 ), having 25 times the 
global warming potential of C0 2 on a 100-year 
timescale (Forster et ah, 2007). Arctic and especially 
high-Arctic regions are already exposed to, and 
predicted to experience, a strong temperature 
increase until the end of the century (4-8 °C higher 
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annual surface air temperatures) (IPCC, 2007), 
which is expected to lead to expanded frost-free 
vegetation periods and increased active layer depths 
in permafrost soils (Tarnocai et al., 2009). 

The major factors in SOC degradation in peat- 
lands, like in mineral soils, are microorganisms of 
the bacterial, archaeal and eukaryotic (fungi) 
domains of life, participating in a cascade of aerobic 
and anaerobic degradation steps, eventually result- 
ing in the emission of CH 4 and C0 2 . Major parts of 
the peat carbon are plant polymers such as cellulose 
and hemicellulose. The degradation of these to 
oligomeric and monomeric sugars is considered 
one of the most important steps in anaerobic 
degradation of SOC (Kotsyurbenko, 2005) and is 
catalysed by a diverse set of hydrolytic extracellular 
enzymes produced by microorganisms. Further 
steps in the anaerobic degradation are carried out 
through anaerobic respiration (like denitrification), 
fermentative and methanogenic pathways, while 
microbial methane oxidation constitutes the biolo- 
gical filter for methane emissions from peat. Phe- 
nolic substances accumulate to high concentrations 
in the anaerobic, water-logged layers of peat soils, 
partly owing to low activity of phenol oxidases, 
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which require oxygen for function. The inhibiting 
effect of phenolic substances has been suggested as a 
major factor for the low SOC degradation rates in 
peat soils (Fenner and Freeman, 2011). 

Several studies have targeted the microorganisms 
involved in CH 4 formation (for example, Ganzert 
et al, 2007; Kotsyurbenko et al, 2007; Metje and 
Frenzel, 2007; H0j et al, 2008) and oxidation (for 
example, Wartiainen et al, 2003; Liebner et al, 
2009; Martineau et al, 2010; Graef et al, 2011) in 
Arctic soils. Broader studies have also been con- 
ducted, for example, on the biogeography of soil 
bacterial communities in different climatic zones 
(Chu et al, 2010), the dynamics of Arctic soil micro- 
bial communities in relation to the composition of 
plant communities (Zak and Kling, 2006) and the 
availability of SOC (Waldrop et al, 2010 and Coolen 
et al, 2011). In recent years, metagenomics and 
metatranscriptomics have developed into powerful 
tools in microbial ecology that enable non-targeted 
studies of the genetic potential, gene expression and 
the composition of soil microbial communities (for 
example, Tringe et al, 2005 and Urich et al, 2008). 
Until now only two metagenomic studies on Arctic 
soils have been conducted, comparing the commu- 
nity in the permafrost with its overlaying active 
layer (Yergeau et al, 2010) and investigating the 
response of a permafrost microbial community to 
thaw (Mackelprang et al, 2011). 

It is in the active layers, where the SOC degrada- 
tion in Arctic peat soils occurs, that more knowledge 
about the processes and organisms involved is 
required to understand these and to predict the 
magnitude of CH 4 and C0 2 emissions from Arctic 
peat in a warmer climate. We have in this study 
applied, for the first time, a combined metagenomic 
and metatranscriptomic approach on high-Arctic 
peat soils to obtain on-site information about the 
genomic potential for SOC transformation and to 
identify the active microorganisms driving these 
processes. We compared the Arctic peat metagenomes 
with metagenomes from temperate soils and proposed 
a model of SOC degradation in Arctic peat, including 
anaerobic respiratory and fermentative pathways, and 
methanogenic and methanotrophic microorganisms. 
The data allow some predictions about future CH 4 
and C0 2 emissions from Arctic peat soils. 

Materials and methods 

Study sites and sampling 

We have investigated the active layers of two high- 
Arctic fens on Svalbard (Norway). Solvatn 
(N78°55.550 Ell°56.611) and Knudsenheia (N78°56. 
544 Ell°49.055) were sampled in August 2009 at 
the peak of the growing season (Supplementary 
Figure Si and S2). Both sites are located in the 
vicinity of the research settlement Ny-Alesund. 
Solvatn is situated on a marine terrace right next 
to the settlement, while Knudsenheia is located 
~5 km northwest of Ny-Alesund. For both sites, 



Knudsenheia and Solvatn, two biological replicates 
were prepared. Three peat blocks (20 x 30 x 20 cm) 
were cut from each of the replicate sites, Si and S2 
(Solvatn) and Kl and K2 (Knudsenheia; see 
Supplementary Figures S1-S3), and transported to 
the laboratory within 2 h after sampling. The active 
layer depths was ~40 cm at both sites. The peat 
blocks were subsequently stored for 2-3 h at 4 °C 
until further processing. The three blocks were 
separated into upper (more oxic) and lower (more 
anoxic) layers, and the corresponding layers, within 
each replicate site, were pooled to account for 
heterogeneity of each replicate site (Supplementary 
Figure S3). The layer separation was based on visual 
distinction of the layers based on its colour and 
structure. The lower layers were processed under 
nitrogen atmosphere to avoid extensive oxygen 
contamination. After pooling, samples were imme- 
diately shock frozen in liquid nitrogen. The whole 
procedure took ~5-10min. Soil pH was measured 
using the KC1 extraction method on fresh peat 
samples. Water content was estimated by drying 
peat over night at 120 °C and reweighing the 
samples. To estimate organic matter contents, the 
dried samples were combusted at 450 °C overnight. 
The samples were reweighed, and the ash weight 
was calculated and added to the weight of organic 
matter. Concentrations of volatile fatty acids and 
ethanol were determined by high pressure liquid 
chromatography analysis of pore water (Metje and 
Frenzel, 2005) extracted from soil samples before 
freezing. Results were integrated using peak simple 
version 393, and ImM standards. 

Nucleic acids extraction, reverse transcription and 454 
sequencing 

From each layer of the biological duplicates 
(Supplementary Figure S3), six or more parallel 
extractions of nucleic acids were performed using a 
modified version of the Griffith's protocol (Urich 
et al, 2008). In samples for metatranscriptomic 
analysis, DNA was digested using the RQl DNase 
treatment (Promega, Madison, WI, USA), followed by 
RNA purification using the MEGAclear kit (Ambion, 
Austin, TX, USA). Samples were subjected to reverse 
transcription using the Superscript II double-strand 
complementary DNA synthesis kit (Invitrogen, Carlsbad, 
CA, USA), following the manufacturer's protocol, 
with the exception that both first- and second-strand 
synthesis was carried out for 4h. RNA template 
addition was in the range of 500-1000 ng. 
Concentrations of double-stranded complementary 
DNA were estimated by SYBR Green I (Invitrogen) 
assay (Leininger et al, 2006). Nucleic acids for 
metagenomic analysis were subjected to RNase A 
treatment for 20min at 37 °C, followed by phenol:- 
chloroform:isoamylalcohol extraction and chloroform: 
isoamylalcohol extraction. DNA concentrations 
were estimated using a NanoDrop spectrophot- 
ometer (Thermo Fisher Scientific, Madison, WI, 
USA). Roche 454 GS FLX Titanium sequencing 
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(454 Life Sciences, Branford, CT, USA) was carried 
out at the CEES at the University of Oslo. 

Bioinformatic analysis 

Metatranscriptomic and metagenomic sequences 
were first filtered using LUCY (Chou and Holmes, 
2001), removing short (<150bp) and low-quality 
sequences (>0.2% error probability) . Metatranscrip- 
tomic sequences stemming from ribosomal RNA 
(rRNA) and putative messenger RNAs (mRNAs) 
were separated by comparing all sequences against 
a combined database of small and large subunits of 
rRNA (SSU and LSU rRNA) using BLASTN 
(Altschul et ah, 1997) and MEGAN (Huson et aL, 
2007; Urich et aL, 2008). Sequences with a bit score 
<50 were assigned as putative mRNA tags. Metage- 
nomic and putative mRNA tags in the metatran- 
scriptomes were subjected to the 454 Replicate filter 
(Gomez-Alvarez et aL, 2009) for the removal of 
artificially replicated sequences. The filter was 
applied to remove exact duplicate sequences from 
the metatranscriptomes (settings: sequence identity 
cutoff: 1.0, length difference requirement: 1.0 and 
number of beginning base pairs to check: 3). For the 
metagenomes, the default settings were applied. 

SSU ribo-tags were taxonomically assigned by 
MEGAN analysis of a BLASTN file against a SSU 
rRNA reference database (parameters: min. bit score 
150, min. support 1, top percent 10; 50 best blast 
hits) (Urich et aL, 2008; Lanzen et aL, 2011). SSU 
ribo-tags assigned to the order Methylococcales were 
assembled into ribo-contigs using CAP3 (Huang and 
Madan, 1999), using two subsequent rounds of 
assembly with (1) a minimum overlap of 150 bp 
with a minimum similarity threshold of 99% and 
mismatch and gap scores of — 130 and 150, and (2) 
minimum overlap of 150 bp and minimum 97% 
similarity threshold, respectively. Contigs were 
uploaded in ARB (http://www.arb-home.de/), and 
contigs longer than 1300 bp together with selected 
reference sequences were used to generate a max- 
imum likelihood tree with the RAxML (Stamatakis 
et aL, 2005) algorithm implemented in ARB using 
default settings. 

The metagenome sequences and the mRNA frac- 
tion of the metatranscriptome were functionally 
annotated for initial screening using the metage- 
nomics (MG)-RAST server with maximum e- value 
for a significant match set to le — 4 (Meyer et aL, 
2008). The sequences were also taxonomically 
binned by MEGAN analysis (parameters: min. bit 
score 50, min. support 1, top percent 2; 50 best blast 
hits) of BLASTX files against the RefSeq protein 
database (e-value <le — 4). Screening for specific 
functional genes and transcripts were carried out 
using custom reference databases generated from the 
UniProtKB/Swiss-Prot or UniProtKB/TrEMBL data- 
bases. BLASTX searches were carried out with an 
e-value threshold of le — 10. All sequences below 
threshold were selected as query sequences for a 
BLASTX search with the same parameters against 



the RefSeq database. Only those that were assigned 
the same function in the second round were 
considered positives. 

In the metagenomes, the genes encoding carbohy- 
drate-active enzymes were annotated. The metage- 
nomic reads were translated into all six frames, each 
frame into separate ORFs (open reading frames), 
avoiding any '*' characters marking stop codons in a 
resulting ORE All ORFs equal to 40 amino acids or 
larger were screened for assignable conserved 
protein domains. Glycoside hydrolases (GH), pro- 
teins with a carbohydrate-binding motif, carbohy- 
drate esterases, polysaccharide lyases and other 
carbohydrate-active enzymes were inspected by 
reference HMMs (Hidden Markov Models) using 
HMMER tools (http://hmmer.janelia.org/) with the 
PFAM database HMMs (PFAM release 25, http:// 
pfam.janelia.org). All database hits with e- values 
below a threshold of 10-4 were counted. For 
the GH44 family, no PFAM HMM is available. 
Therefore, representative sequences were selected 
from the CAZy website (http://www.cazy.org), the 
sequence regions corresponding to the family were 
determined and used to create a HMM for the GH44 
family. The screening was performed on the HPC 
computer STALLO at the University of Troms0 
(http://docs.notur.no/uit). A selection of sequences 
was also used as query sequences in BLASTP 
searches against the RefSeq database of the NCBI. 
Within the chosen threshold of the HMMER, all 
identified carbohydrate-active enzymes gave either a 
corresponding hit or no hit in the RefSeq database. 
For the taxonomic assignment of gene and transcript 
sequences, the corresponding BLAST outputs were 
uploaded in MEGAN (parameters: min. bit score 35, 
min. support 1, top percent 10; 50 best blast hits). 

Statistical analysis 

Significant differences between the frequencies of 
conserved protein domains in ORFs of different soils 
were evaluated statistically by using the R package 
(R_Development_Core_Team, 2009) using the y 2 con ~ 
tingency table test. The contingency table holds the 
frequency counts of hits and non-hits for a certain 
PFAM domain category of two different soils. The 
total frequency count is given by all hits found for any 
domain in the PFAM database. In cases where the 
frequencies are too low to meet the rules of the test, 
the probabilities of observing a sample statistic as high 
as the test statistic were calculated by Monte Carlo 
simulations with 100 000 replicates. 

Data deposition 

The sequence data generated in this study was 
deposited in the Sequence Read Archive of NCBI and 
are accessible through accession number SRP014474. 

Results and discussion 

The active layers of the two permafrost fens 
at Solvatn and Knudsenheia were slightly acidic 
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(pH5-6), had a water content of -70-90%, an 
organic matter content of 40-90% (Supplementary 
Table Si) and were characterised by a moss cover 
dominated by Calliergon richardsonii (Solheim 
et al, 1996; Wartiainen et al, 2003; H0j et al, 
2005). The mosses were interspersed by grasses 
[Dupontia pelligera), which were heavily sup- 
pressed by grazing Barnacle geese. High concentra- 
tions of acetate, lactate and ethanol indicated 
anaerobic metabolisms through fermentative and 
acetogenic pathways, even in the top layers 
(Supplementary Table Si). Remarkably high yields 
of nucleic acids (DNA and RNA) per gram dry peat 
were obtained from all samples and layers 
(Supplementary Table S2), which suggests a high 
microbial abundance and activity at the peak of the 
growing season. The yields, especially of RNA, 
decreased with depth, indicative of decreasing 
microbial activity. Replicated metatrans crip tome 
data sets were obtained from all upper layers 
(Supplementary Figure S3) whereas the reverse 
transcription was severely inhibited in the lower 
layers, presumably owing to phenolic substances 
(Supplementary Table S2). Metagenomic libraries 
were obtained from all layers of one duplicate from 
Solvatn and Knudsenheia (Kl and S2, Supplement- 
ary Figure S3). 454 Titanium pyrosequencing of the 
metatranscriptomes and metagenomes resulted in 
120 000-190 000 sequences per library, with one 
exception (S2b complementary DNA, Supplement- 
ary Table S3). The number of rRNA sequences in 
these libraries ranged from 100 000 to 150 000 while 
the number of putative mRNA sequences ranged 
from 7000 to 8500. 



Community structure of active microorganisms 
The community composition of bacteria, archaea 
and eukaryotes, as determined from the small 
subunit rRNA fragments (SSU ribo-tags) of the 
metatranscriptomes (up to 74 000 ribo-tags), 
revealed that the biological replicates of both 
sites had a similar composition of taxa at 
domain, phylum and class levels of resolution 
(Supplementary Figure S4). Ribo-tags are consid- 
ered a measure of living biomass, as they stem from 
ribosomes, most likely indicative of active organ- 
isms (Urich et al, 2008; Urich and Schleper, 2011). 
Bacteria comprised most ribo-tags; — 70%-80% in 
all upper layers and —84% in the lower, anoxic 
layer of Solvatn, S2b (Figure 1). Eukaryotic ribo-tags 
comprised ~ 15%— 30%, whereas the fraction of 
Archaea was surprisingly small in all the upper 
layers (~0.01%-0.13%), but higher in S2b (-1%). 
The majority of archaeal sequences stemmed from 
methanogenic archaea. Similar low relative abun- 
dances of Archaea have been reported from Alaskan 
permafrost soils (for example, Waldrop et al., 2010; 
Wilhelm et al, 2011). The bacterial community 
composition in both peat soils was very similar, 
with most ribo-tags stemming from eight phyla, most 



notably the Proteobacteria (~ 37%-45% of bacterial 
ribo-tags), Actinobacteria ( — 15%), Planctomycetes 
(~9%-14%), Verrucomicrobia (~9%-ll%), Acid- 
obacteria (~8%-10%) and the Chloroflexi (~7%- 
9%). Within the Proteobacteria, Deltaproteobacteria 
(mainly Myxococcales) were the most abundant 
class (~19%-21% of bacterial ribo-tags). Alphapro- 
teobacteria made up —8%, Betaproteobacteria 
comprised —5% while Gammaproteobacteria repre- 
sented on average 3.2% of the total community. The 
major eukaryotic taxon was the Protist kingdom 
Alveolata, putative grazers of bacteria, which repre- 
sented ~17%-20% of eukaryotic ribo-tags in both 
layers of Solvatn and 6% -8% in the top layer of 
Knudsenheia (Figure 1). Other abundant Protists 
were the Amoebozoa, Rhizaria (both predators of 
bacteria) and Stramenopiles. Interestingly, the Pro- 
tists constituted the major eukaryotic group, exceed- 
ing even the fraction of fungal (~2%-8%) and 
metazoan (~11%-13%) ribo-tags (Figure 1). This is 
in contrast to temperate and boreal mineral soils, 
where the fungi and metazoa comprise the majority 
of either eukaryotic ribo-tags (Urich et al., 2008) or 
biomass (Schaefer, 1990; Berg and Bengtsson, 2007) 
and the Protists are much less abundant. Our data 
indicate that the bacterial energy channel might be 
more important than the fungal energy channel in 
these peat soils. Metazoan ribo-tags originated from 
typical micro- and meso-faunal groups of higher 
trophic levels, including Nematoda, Tardigrada, 
Rotifera and Arthropoda. Ascomycota, and to a 
minor extent Basidiomycota, dominated the small 
fungal community. At Solvatn, a change in abun- 
dance with depth was observed for several taxa 
(Figure 1). While Alphaproteobacteria, Betaproteo- 
bacteria and Chloroflexi were present in larger 
fractions in the lower layer, Planctomycetes and 
Verrucomicrobia decreased with depth. Most notably, 
however, was the increase in the fraction of Archaea 
and the decrease of Protists and fungi with depth. The 
ribo-tags originated mostly from actively transcribing 
organisms, indicated by a taxonomic binning of the 
mRNA fraction, which showed similar community 
profiles for the bacterial taxa (Figure 2a). 

There are few studies addressing Arctic peat soils 
to compare our metatranscriptomic data with. 16S 
rRNA gene clone library, and metagenomic studies 
of permafrost mineral soils have shown that the 
Actinobacteria, Proteobacteria and Acidobacteria 
dominate these environments (Hansen et al., 2007; 
Liebner et al, 2008; Steven et al, 2008; Chu et al, 
2010; Yergeau et al, 2010; Wilhelm et al, 2011), 
while clone library, T-RFLP and FISH studies of a 
sub-Arctic acidic sphagnum peat report a domi- 
nance of Alphaproteobacteria, Betaproteobacteria, 
Acidobacteria, Verrucomicrobia and Planctomycetes 
(Dedysh et al, 2006; Pankratov et al, 2011). A recent 
study of permafrost peat showed that the active 
layer communities were dominated by Actino- 
bacteria, Proteobacteria and Chloroflexi, while 
the permafrost microbiota also contained large 
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Figure 1 Three-domain community profile of the microbiota in Svalbard peatlands. The figures are based on the ribo-tags fraction of the 
metatranscriptome. The size of the boxes is proportional to the fraction of ribo-tags of the respective taxa. (a) Top layer of Solvatn peat 
generated from two biological replicate data sets (Sla and S2a). (b) Lower layer of Solvatn generated from one data set (S2b). (c) Top layer 
of Knudsenheia peat generated from two biological replicate data sets (Kla and K2a). Asco, Ascomycota; Basidio, Basidiomycota; 
Meth.bac, Methanobacteriales', Meth.mic, Methanomicrobiales; Meth.sarc, Methanosarcinales. 



populations of Bacteroidetes and Firmicutes 
(Mackelprang et al., 2011). Interestingly, taxa such as 
the Verrucomicrobia and Planctomycetes, which were 
among the most abundant taxa in the Svalbard peat 
soils as well as the acidic Sphagnum peat, were not 
reported as major taxa in the studies on mineral soils. 



Genomic potential for degradation of plant polymers 
Plant polymer and phenolic compound degradation 
are key processes of SOC decomposition in peat 
soils. We analysed the genomic potential for these 
processes, by screening the metagenomes for genes 
encoding conserved protein family (PFAM) domains 
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Figure 2 Taxonomic assignment of metatranscriptomic and metagenomic sequences. The community structure is displayed at phylum 
level resolution with Proteobacteria split into classes. Sequences assigned to Proteobacteria refer to sequences that could not be assigned 
to class level resolution. S2 and Kl indicates where the samples were collected, in Solvatn and Knudsenheia respectively, while a, b and 
c indicates the depth of the sample from top and down (Supplementary Figure S3), (a) Taxonomic assignment of ribo-tags, mRNA and 
metagenomic DNA (gDNA) to the domain bacteria, (b) Taxonomic assignment of the metagenomic sequences encoding polysaccharide- 
degrading enzymes (cellulases, endohemicellulases and debranching enzymes). All sequences assigned to these categories were pooled 
together and taxonomically binned using MEGAN (see Materials and methods for details). 



of relevant enzymes, that catalyse the hydrolysis of 
abundant plant polymers such as cellulose and 
hemicelluloses, and the degradation of phenolic 
compounds (for example, lignin) (Table 1). Despite 
the rather uniform plant cover, dominated by the 
moss C. richardsonii, we identified a large variety of 
CAZy (carbohydrate-active enzyme) families, 
including 76 GH families, 36 carbohydrate-binding 
module, 3 phenol oxidase, 3 polysaccharide lyase 
and 5 carbohydrate esterase families (Supplement- 
ary Table S4). GH5 were the most abundant among 
the families of cellulases. Genes encoding endohe- 
micellulases (especially GH families 10 and 26 
that target xylan), debranching enzymes and 
phenol oxidases involved in the degradation of 
phenolic compounds were also abundant (Table 1). 

We wanted to know if the microbiotas of our 
moss-dominated, carbon-accumulating high-Arctic 
peat soils contain a different genetic potential for the 
degradation of plant polymers and recalcitrant 
phenolic compounds than the microbiota of soils 
from other climate zones, having a vascular plant 
cover. Therefore, we compared our metagenomic 
data sets with metagenomes from temperate grass- 
lands and farmlands, subtropic rain forest soils and 
Arctic permafrost mineral soils (see Supplementary 
Table S5 for details about the soil data sets), 
focussing on the PFAMs relevant for plant polymers 
and phenolic compound degradation. The functional 
potential was similar in all soils, and genes for the 
same PFAMs dominated the metagenomes (Supple- 
mentary Table S6). A x 2 test did not detect statistically 
significant differences between the metagenomes from 
arctic peat and the other metagenomes (Supplement- 
ary Table S7), regardless of the variations in geography 
or plant coverage. 

The chemical composition of moss cell walls is 
generally considered to be simpler than the one from 



vacular plants; mosses have a different and less 
complex hemicellulose composition, and do not 
contain lignin, but a simpler and less abundant 
phenolic compound called lignan (Pena et al., 2008; 
Sarkar et al., 2009; Popper et al., 2011), with the 
exception of Sphagnum species containing an 
abundant lignin-like polymer (Ligrone et al., 2008). 
It was therefore surprising that the microbiotas 
harboured such a similar genomic potential. This 
result might to some extent be explained by 
methodological constraints. Fon instance, the 
PFAMs for phenolic compound degradation do not 
distinguish between lignin and lignan as substrate. 
Also, the enzymes themselves might have a broad 
substrate specificity. Finally, the presence of many 
similar types of polymers in both mosses and 
vascular plants should be reflected in an overall 
similar genetic repertoire for degradation in the 
respective soil microbiotas. It might well be that the 
actual composition of plant polymers is only 
reflected in the gene expression patterns of the 
microbiota and not in its genetic potential. However, 
the rather small number of mRNAs in the metatran- 
scriptome prevented us from studying this. 

Another major constraint for SOC degradation, 
besides the short growing season, is the accumula- 
tion of phenolic compounds in the water-logged, 
predominantly anaerobic peat soils, which inhibits 
enzymatic activity (Fenner and Freeman, 2011). 
Expression of phenol oxidases was detected in the 
top layer of Solvatn. However, the genetic potential 
for phenol oxidases was observed in all peat layers, 
indicating that a degradation of these inhibitory 
compounds can occur if the peat soils get oxyge- 
nated, possibly through more pronounced drought- 
rewetting cycles at higher temperatures, like it has 
recently been shown for temperate peatlands 
(Fenner and Freeman, 2011). 
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Table 1 PFAMs that target plant polymers, identified in metagenomes from high-Arctic peatlands of Svalbard 





oza 




Kia 


KID 


KlC 


V_jC 11 LllLloC O 












GH5 


4 (2.3) 


4 (1.8) 


4 (2.1) 


8 (3.9) 


18 (5.8) 


GH6 


2 (1.2) 


0 (0) 


1 (0.5) 


1 (0.5) 


2 (0.7) 


GH7 


0 (0) 


1 (0.5) 


0(0) 


1 (0.5) 


0 (0) 


GH9 


1 (0.5) 


2 (0.9) 


4 (2.1) 


1 (0.5) 


6 (2.0) 


GH44 


2 (1.2) 


2 (0.9) 


3 (1.6) 


2 (1.0) 


0(0) 


GH45 


1 (0.6) 


0(0) 


0(0) 


0(0) 


1 (0.3) 


GH48 


0 (0) 


0(0) 


0(0) 


0(0) 


0(0) 


Sum 


10 (5.9) 


Q f4 1 1 


12 (6.3) 


13 (6.3) 


27 (8.8) 


Debranching enzymes 












GH51 


5 (2.9) 


8 (3.6) 


3 (1.6) 


5 (2.4) 


13 (4.2) 


GH54 


1 (0.6) 


3 (1.4) 


0(0) 


0(0) 


0(0) 


GH62 


0 (0) 


0(0) 


1 (0.5) 


0(0) 


0(0) 


GH67 


2 (1.2) 


1 (0.5) 


0(0) 


0(0) 


2 (0.7) 


GH78 


1 (0.6) 


2 (0.9) 


2 (1.1) 


10 (4.9) 


15 (4.9) 


Sum 


9 (5.3) 


14 (6.4) 


6 (3.1) 


15 (7.3) 


30 (9.7) 


J^nrinhpmirplln ln<zp<z 

1—1 1 L vl KJ1 L I'll 1 1 O K ' 1 1 UiUiJUO 












GH8 


0 (0) 


0 (0.0) 


1 (0.5) 


0(0) 


2 (0.7) 


GH10 


2 (1.2) 


6 (2.7) 


1 (0.5) 


5 (2.4) 


7 (2.3) 


GH11 


0 (0.0) 


0 (0.0) 


0 (0.0) 


0(0) 


1 (0.3) 


GH12 


0 (0.0) 


1 (0.5) 


0 (0.0) 


0(0) 


0(0) 


GH26 


4 (2.3) 


4 (1.8) 


0(0) 


1 (0.5) 


12 (3.9) 


GH28 


0 (0.0) 


1 (0.5) 


4 (2.1) 


2 (1.0) 


8 (2.6) 


GH53 


0 (0.0) 


3 (1.4) 


3 (1.6) 


2 (1.0) 


4 (1.3) 


Sum 


6 (3.5) 


15 (6.8) 


9 (4.7) 


10 (4.9) 


34 (11.0) 


Oligosacch arid e-d egrading en zym es 












GH1 


9 (5.3) 


10 (4.6) 


6 (3.1) 


14 (6.8) 


20 (6.5) 


GH2 


5 (2.9) 


0 (0.0) 


4 (2.1) 


4 (1.9) 


15 (4.9) 


GH3 


5 (2.9) 


Q f4 1 1 




Q f4 4l 


18 (5.8) 


GH29 


5 (2.9) 


4 (1.8) 


4 (2.1) 


11 (5.3) 


13 (4.2) 


GH35 


2 (1.2) 


0 (0.0) 


0 (0.0) 


0 (0.0) 


3 (1.0) 


GH38 


0 (0.0) 


0 (0.0) 


1 (0.5) 


3 (1.5) 


4 (1.3) 


GH39 


2 (1.2) 


5 (2.3) 


1 (0.5) 


2 (1.0) 


7 (2.3) 


GH42 


3 (1.8) 


2 (0.9) 


2 (1.1) 


6 (2.9) 


11 (3.6) 


GH43 


8 (4.7) 


7 (3.2) 


5 (2.6) 


8 (3.9) 


10 (3.2) 


GH52 


0 (0.0) 


0 (0.0) 


0 (0.0) 


2 (1.0) 


0 (0.0) 


Sum 


39 (23) 


37 (17) 


29 (15) 


59 (29) 


101 (33) 


Phenolic compound-degrading enzymes 












Laccase 


1 (0.6) 


8 (3.6) 


7 (3.7) 


11 (5.3) 


12 (3.9) 


Dioxygenase 


2 (1.2) 


8 (3.6) 


3 (1.6) 


2 (1.0) 


5 (1.6) 


Peroxidase 


4 (2.3) 


5 (2.3) 


3 (1.6) 


8 (3.9) 


1 (0.3) 


Sum 


7 (4.1) 


21 (9.6) 


13 (6.8) 


21 (10.2) 


18 (5.8) 


Tot. ass. seq. 


17067 


21980 


19104 


20581 


30 834 



The metagenomic profile is divided into five categories according to the functional role of the protein families in plant polymer degradation, 
a format previously applied in Allgaier et al. (2010) and Pope et al. (2010). Numbers reflect the number of significant hits (see Materials and 
methods for details). Shown in brackets are the numbers in °/ 00 o of the total number of sequences assigned to a PFAM. Tot. ass. seq., the number of 
sequences that were assigned to a PFAM (See Materials and methods for details). See Supplementary Table S4 for PFAM accession numbers and 
description of protein families. 
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We aimed to identify the microorganisms posses- 
ing the genes encoding these hydrolytic enzymes by 
using taxonomic binning with the MEGAN software 
(See Materials and methods for details). However, 
this type of taxonomic binning is problematic, 
owing to the unequal representation of reference 
genomes for many taxa in the public databases 
(Urich et al., 2008) and horizontal gene transfer. We 
therefore compared the community profiles from 
SSU ribo-tags with the taxonomically assigned 
metagenomic and metatranscriptomic sequences 
from the peat top layers using MEGAN (Figure 2a). 



These profiles were similar for many taxa, suggest- 
ing that genes and transcripts can be used as 
taxonomic markers with some level of confidence. 
The genes assigned to the three major categories of 
polysaccharide degradation (cellulases, endohemi- 
cellulases and debranching enzymes; see Table 1) 
were taxonomically annotated. The bacterial phyla 
Bacteroidetes, Actinobacteria and Verrucomicrobia 
possessed the majority (>70%) of these genes 
(Figure 2b). The profiles of the two peat soils 
differed, with the Knudsenheia metagenomes hav- 
ing a higher fraction of genes from Bacteroidetes and 
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Table 2 Metagenomic and metatranscriptomic reads assigned to enzyme classes that catalyse key steps in anaerobic respiration and 
fermentation, and fermentative pathways 



Pathway/ function E.G. Gene Metagenomes Metatranscriptom.es 



S2a S2b Kla Klb Klc Sla + S2a Kla + K2a 



Nitrate reduction 


1.9.6.1 


nap A 


1 


Nitrate reduction 


1.7.99.4 


narG 


10 


Nitrite reduction 


1.7.2.1 


nirK 


1 


Nitrite reduction 


1.7.2.1 


nirS 


0 


Nitric oxide reduction 


1.7.99.7 


norB 


1 


Nitrous oxide reduction 


1.7.99.6 


nosZ 


3 


Sulphate reduction 


1.8.99.3 


dsrA 


0 


Fermentation: [FeFe] hydrogenase 


1.12.7.2 


hydA 


2 


Homoacetogenesis 


6.3.4.3 


fthfs 


6 


Methanogenesis 


2.8.4.1 


mcrA 


0 


Acetoin, butanediol metabolism 






55 


Acetone, butanol, ethanol synthesis 






243 


Acetyl-CoA fermentation to butyrate 






310 


Butanol biosynthesis 






235 


Fermentations: lactate 






52 


Fermentations: mixed acid 






107 


Tot. ass. sequences 






69176 



6 


7 


4 


1 


o 


o 


5 


14 


19 


28 


o 


o 


1 


2 


3 


3 


0 


0 


1 


0 


2 


1 


0 


2 


4 


3 


2 


4 


0 


1 


3 


1 


2 


1 


0 


0 


1 


1 


4 


11 


0 


0 


9 


0 


7 


8 


0 


0 


12 


7 


10 


23 


0 


0 


0 


1 


0 


2 


0 


0 


92 


60 


85 


133 


8 


8 


298 


319 


280 


346 


17 


17 


362 


406 


329 


407 


16 


22 


293 


315 


263 


308 


9 


19 


76 


56 


74 


104 


2 


3 


157 


121 


131 


190 


5 


8 


63 527 


85 664 


54 889 


72 166 


6315 


6205 



First column shows function or main pathway, second column shows the enzyme class entry of the ExPASY bioinformatics resource, while third 
column shows the gene abbreviation. Merged biological replicates for the two sites Solvatn (Sla + S2a) and Knudsenheia (Kla + K2a) are listed 
under metatranscriptomes. Total ass. sequences, the number of sequences assigned a function by MEGAN analysis. Sequences assigned to 
pathways were annotated using (MG)-RAST (see Materials and methods for details). 



Actinobacteria, while Solvatn had a higher fraction 
of genes assigned to Verrucomicrobia. Surprisingly, 
few genes were taxonomically assigned to fungi 
(Figure 2b). This was in accordance with the low 
abundance of fungal, and especially basidiomycotal 
ribo-tags (Figure 1). The latter have been reported as 
particularly potent decomposers of SOC (Waldrop 
et al, 2010). Low fungal abundances in permafrost 
soils and the influence of plant community compo- 
sition on fungal population size have been shown 
before (Zak and Kling, 2006; Yergeau et al, 2010; 
Waldrop et al, 2010). 

Anaerobic respiration and fermentation 
We investigated the genomic potential for and the 
expression of key enzymes in anaerobic respiration 
and fermentation, processes which drive the anae- 
robic decomposition of SOC. Especially the genes 
encoding dissimilatory nitrate reductases were 
abundant (Table 2). This suggests that denitrifica- 
tion is important, which is supported by the high 
nitrate concentrations detected (70 and 84|ig per 
gram soil in Solvatn and Knudsenheia, respectively; 
Alves, 2011). An important role of denitrification 
has also been suggested by recent metagenomic 
studies of permafrost affected soils (Yergeau et al, 
2010; Mackelprang et al, 2011). Despite the high 
abundance of denitrification genes, only nirS (nitrite 
reductase) and norB (nitric oxide reductase) tran- 
scripts were identified (Table 2); however, the small 
size of the mRNA pools in the metatranscriptomic 
data sets prevented further statements. Also sulphate 
reduction might have a role, as genes of the key 
enzyme DSR were found and their abundance 



increased with soil depth in Knudsenheia (Table 2). 
The majority of the identified anaerobic respiration 
genes were assigned to taxa within Actinobacteria and 
Proteobacteria (Supplementary Figure S5), although 
the genes were broadly distributed among several 
different taxa, as previously described for denitrifyers 
(Philippot and Hallin, 2005) and sulphate reducers 
(Wagner et al, 2005). Genes and transcripts of several 
fermentative pathways were detected (Table 2), 
explaining the high concentration of fermentation 
products in the peat soils (Supplementary Table Si). 
Likewise, genes encoding the terminal hydrogenase of 
H 2 -evolving fermentations [hydA) were found. Their 
abundance increased with depth, as did the abun- 
dance of formyltetrahydrofolate synthetase genes [fhs], 
which encodes the key enzyme of the Acetyl-CoA 
pathway of homoacetogenesis. The majority of hydA 
and fhs genes were assigned to Firmicutes and 
Actinobacteria, suggesting an important role of taxa 
within these phyla in the anaerobic degradation of 
SOC in the Svalbard peat (Supplementary Figure S6). 
Taken together, these findings indicate a gradual 
transition from aerobic to anaerobic degradation of 
SOC with depth. 



Methane cycling 

The fraction of methanogenic archaea increased with 
depth (Figure 1 and Supplementary Figure S6), 
reflecting the presumably more anoxic conditions 
and higher concentrations of fermentation inter- 
mediates (Supplementary Table Si). Three orders of 
methanogenic archaea were found, Methanobacter- 
iales, Methanomicrobiales and Methanosarcinales, 
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narl 

I 



narH 

I 



Nitrite 
transporter 




/ 

norJ 



narG 



Hyp. Prot 
nirC \ nirDnirG nirJ 

DimmE> 

/ l I I I 
nirS nirM nirF nirL nirH 



■ Contig 21, Knudsenheia, 1358, 16 

I r Contig 1 0, Knudsenheia, 1 343, 28 

L Contig 3, Solvatn, 1 377, 40 
j- Contig 7, Knudsenheia, 1435, 438 
1— Contig 4, Solvatn, 1 350, 28 
Contig 2, Solvatn, 1366, 14 



■r 



Contig 4, Knudsenheia, 1450, 309 
Contig 9, Knudsenheia, 1 330, 81 



Contig 14, Knudsenheia, 1430, 176 
— Contig 2, Knudsenheia, 1 330, 59 
Lr Methylobacter psychrophilus, 1527, AF 152597 
n Methylobacter tundripaludum, 1487, AJ41 4655 

j — Contig 18, Knudsenheia, 1342, 37 
J 1 Contig 3, Knudsenheia, 1443, 391 

Contig 1 , Solvatn, 1 382, 28 
p Contig 12, Knudsenheia, 1418,89 

L Contig 19, Knudsenheia, 1361,88 

-j Crenothirx polyspora 
- Methylosoma difficile, 1491, DQ1 79050 




- Contig 20, 98, Knudsenheia, 1 382, 36 
Methylobacter marinus, 7 499, AF304 7 97 
— Methylobacter whittenburyi, 1356,120843 
Methylobacter luteus, 1500,AF304195 

Methylomicrobium 
Methylosarcina 



Methylocycstis parvus, 1408, Y 18945 
Methylocapsa acidiphila, 1415, AW 8726 
Methylocella silvestris, 1409, AJ491847 
-\J Nitrobacter 



Methylomonas 
Methylosphaera hansonii, 1451, U67929 
C7 Methylocaldum 

Methylococcus 

Methylohalobius crimeensis, 1 462, AJ58 1 837 

Methylothermus thermalis, 1 462, AY829009 



0.10 



Figure 3 M. tundripaludum in Svalbard peat soils. The phylogenetic tree shows assembled SSU rRNA contigs of type I methanotrophs. 
Most of the nearly full-length ribo-contigs (14 out of 16) are >97% identical to M. tundripaludum. Ribo-contig description includes the 
following: contig id, site (Knudsenheia-green, Solvatn-red), contig length and the number of single ribo-tags that went into the assembly. 
The reference sequence description (black) includes sequence length and accession number. The length of the bar indicates 0.10 changes 
per nucleotide. The tree was constructed using the ARB software (See Materials and methods for details). The insert shows the nitrate 
reductase (open reading frames (ORFs) 3936-3940) and nitrite reductase operons (ORFs 3936-3940) identified in the genome of 
M. tundripaludum. 



Table 3 mRNAs assigned to M. tundripaludum in the Knudsenheia metatranscriptome 



Annotation AC Function No. of reads 



pmoC 


ZP 


,07656564. 1 


Methane oxidation 


8 


pmoA 


ZP 


_07656563.1 


Methane oxidation 


7 


pmoB 


ZP 


_07656562.1 


Methane oxidation 


6 


CsbD family protein 


ZP 


_07656622.1 


Stress response 


2 


MgtE intracellular region 


ZP 


_07654467.1 


Magnesium transport 


2 


Transport-associated protein 


ZP 


_07656600.1 


Transport 


2 


Methanol dehydrogenase 


ZP 


_07652313.1 


Methanol oxidation 


1 


6-Phospho 3-hex. 


ZP 


_07653633.1 


RuMP pathway step 2 


1 


Succinate dehydrogenase 


ZP 


_07652392.1 


Electron transport Complex II 


1 


NADH dehydrogenase 


ZP 


_07653777.1 


Electron transport Complex I 


1 


NADH:ubiquinone oxidoreductase 


ZP 


_07655325.1 


Electron transport Complex I 


1 


Elec. trans, prot. SenC 


ZP 


_07652162.1 


Biogenesis of respiratory systems 


1 


Nitrogenase iron protein 


ZP 


_07654124.1 


Nitrogen fixation 


1 



Abbreviations: elec. trans, prot. SenC, electron transport protein SenC; mRNA, messenger RNA; RuMP, ribulose monophosphate. 
Functionally annotated transcripts related to energy, carbon and nitrogen metabolism are shown. mRNA sequences were taxonomically assigned 
to M. tundripaludum by MEGAN (see Materials and methods for details). 
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Figure 4 Schematic overview of the main degradation pathways of plant polymers in the high-Arctic peatlands of Svalbard. The 
pathways are divided into three categories; aerobic (beige), anaerobic (light brown) and processes occurring under both conditions (dark 
brown). Key microbial taxa for the different degradation steps are presented (orange boxes). The figure is adapted from Figure 1 in 
(Conrad, 1999). 



of which the first two are hydrogenotrophic, produ- 
cing CH 4 from H 2 and C0 2 , whereas the Methano- 
sarcinales are metabolically more versatile, carrying 
out hydrogenotrophic, acetoclastic and methylo- 
trophic methanogenesis. These findings are similar 
to previous studies of methanogens from Solvatn 
(H0j et al, 2005; H0j et al, 2006). In Solvatn, the 
increased abundance of acetotrophic methanogens 
with depth correlated with a decrease in the 
concentration of acetate, indicating that the aceto- 
trophic Methanosarcinales acted as the sink for 
acetate. In general, the high concentrations of 
fermentation products like ethanol, especially in 
the lower layers of Knudsenheia, indicated a low 
efficiency in the terminal fermentation steps, possi- 
bly owing to low activity of, or a lack of, established 
interactions between fermentative syntrophic bac- 
teria and their methanogenic counterparts. 

Type I methanotrophs of the order Methylococ- 
cales were abundant in both peat soils (0.1%-2.4% 
of ribo-tags). We did not detect active type II 
methanotrophs, anaerobic methanotrophic archaea 
(ANME), the recently described anaerobic methano- 
troph Candidatus Methylomirabilis oxyfera of the 



NC10 candidate phylum (Ettwig et al., 2010) or the 
Verrucomicrobium Methylacidiphilium (Islam et al., 
2008). We assembled long, sometimes full-length, 
SSU rRNAs from the ribo-tags of the type I 
methanotrophs to determine the phylogenetic posi- 
tion (See Materials and methods (Urich et al., 2008; 
Radax et al., 2012)). The vast majority of the 
assembled ribo-contigs was closely related to 
M. tundripaludum (Figure 3), an arctic methano- 
troph originally isolated from Solvatn (Wartiainen 
et al., 2006). We could taxonomically assign 104 
mRNA tags to genes encoded in the genome of 
M. tundripaludum (Svenning et al., 2011). Although 
this is a small number, it allowed a superficial 
view into the gene expression and metabolism of 
M. tundripaludum in situ. Remarkably, as much as 
33% of the transcripts encoded subunits of the key 
enzyme particulate methane monooxygenase 
(Table 3). This highlights the importance of particu- 
late methane monooxygenase for M. tundripaludum 
cells, as this enzyme catalyses the initial step of both 
the carbon and energy metabolism through the 
oxidation of methane to methanol. Transcripts of 
the methanol dehydrogenase that catalyses the 
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second step of complete oxidation of methane to 
C0 2 and the ribulose monophosphate pathway for 
carbon assimilation, as well as transcripts encoding 
enzymes of the aerobic respiratory chain, were also 
identified (Table 3). Recently, M. tundripaludum 
was shown to be the dominating methanotroph in a 
SIP-RNA experiment with Solvatn peat (Graef et al., 
2011). It was also shown to be the dominant 
methanotroph in two soils from the Siberian and 
the Canadian Arctic (Liebner et al., 2009; Martineau 
et al., 2010). Although these are few studies, they 
point to the importance of methanotrophs closely 
related to M. tundripaludum having a key role in the 
biological CH 4 filter in Arctic soils. The seemingly 
ecological advantage of M. tundripaludum in these 
often water-logged Arctic soils might be owing to its 
N 2 -fixation capability (Wartiainen et al., 2006; 
Svenning et al., 2011), as well as to the presence of 
nitrate reductase [nar] and nitrite reductase [nir] 
operons in its genome (Figure 3) (Svenning et al., 
2011). It is tempting to speculate that M. tundripa- 
ludum could employ the NAR and NIR enzymes in 
denitrifying anaerobic methane oxidation, utilising 
a pathway similar to that of Candidatus M. oxyfera 
(Ettwig et al., 2010). Alternatively, it could have the 
ability to grow anaerobically on simple organic 
compounds using denitrification. Experiments are 
underway to test these hypotheses. 
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for anaerobic bacteria and methanogenic archaea, in 
the cases where the peat soils remain water- 
saturated. Here, CH 4 production likely will increase, 
although this might well be balanced by methano- 
trophic bacteria, of which M. tundripaludum seems 
to be particularly important. However, to which 
extent such key populations can adapt to the 
changing conditions remains a crucial question. 
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Conclusion 

We have described the metabolic potential and 
activity of microorganisms involved in SOC degra- 
dation in high- Arctic peat soils (Figure 4). This first 
combined metagenomic and metatranscriptomic 
study on any soil (to our knowledge), revealed taxa 
that are important to specific processes such as 
hydrolysis of plant polymers (for example, Verruco- 
microbia and Bacteroidetes), fermentations (for 
example, Firmicutes), methanogenesis and metha- 
notrophy. Actinobacteria seemed to be particularly 
important, having a metabolic potential for carrying 
out several of the key steps in SOC degradation 
(Figure 4). Based on their relative abundance, fungi 
appeared to not have an important role. However, as 
warming will lead to an extended growing season, 
this might result in a change of the plant covers from 
mosses to vascular plants and in more oxygenated 
active layers, where fungi might become more 
important for SOC decomposition. Dryer and more 
oxygenated active layers will then likely lead to 
increased peat decomposition, owing to lowered 
concentrations of phenolic substances, caused by 
increased activity of aerobic microorganisms that 
synthesise phenol oxidases (Fenner and Freeman, 
2011). These combined effects might increase the 
rate of SOC degradation above that of carbon 
sequestration, and turn these peat soils into C0 2 
sources. Increased active layer depth and higher soil 
temperatures will likely result in bigger habitat size 
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